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£1 ! Abstract 

o 

The representation of the seesaw orthogonal matrix in the complex plane estab- 
lishes a graphical correspondence between neutrino mass models and geometrical 
Ph. configurations, particularly useful to study relevant aspects of leptogenesis. We 

Ph. first derive the CP asymmetry bound for hierarchical heavy neutrinos and then an 

expression for the effective leptogenesis phase, determining the conditions for max- 
imal phase and placing a lower bound on the phase suppression for generic models. 
^ . Reconsidering the lower bounds on the lightest right-handed (RH) neutrino mass 

Mi and on the reheating temperature T reri , we find that models where one of the 
two heavier neutrino masses is dominated by the lightest right-handed (RH) neu- 
trinos, typically arising from connections with quark masses, undergo both phase 
suppression and strong wash-out such that Mi(T re h) > 10 11 (10 10 )GeV. The win- 
dow 10 9 GeV < Mi, T re h ^ 10 10 GeV is accessible only for a class of models where 
mi is dominated by the lightest RH neutrino, with no straightforward connections 
with quark masses. Within this class we describe a new scenario of thermal lepto- 
genesis where the baryon asymmetry of the Universe is generated by the decays of 
the second lightest RH neutrino, such that the lower bound on Mi disappears and 
is replaced by a lower bound on Mi- Interestingly, the final asymmetry is indepen- 
dent on the initial conditions. We also discuss the validity of the approximation of 
hierarchical heavy neutrinos in a simple analytical way. 



1 Introduction 



The seesaw mechanism [1] is an elegant theoretical explanation for the lightness of neu- 
trinos compared to the other Standard Model particles. Neutrino masses are predicted to 
be small but non-vanishing and this has been successfully confirmed in neutrino mixing 
experiments [2]. Moreover, within the seesaw, neutrino mixing data point to a new scale 
~ M% w I \yj 'Am s 2 ol - y/Amltn) ~ 10 14 ~ 15 GeV, compatible with GUT's expectations [3]. 

Beyond neutrino mixing and masses, the seesaw mechanism has different phenomeno- 
logical implications with interesting testable predictions. An intriguing cosmological con- 
sequence of the seesaw mechanism is leptogenesis [4]. This provides an attractive expla- 
nation for the observed baryon asymmetry of the Universe. Analogously to the case of 
GUT baryogenesis models [5], the asymmetry is generated from heavy particle decays. 
However, in the case of leptogenesis, a B-L asymmetry is generated in the form of lepton 
number that is partly converted into a baryon asymmetry by (B-L conserving) sphaleron 
processes [6] . In a thermal scenario the initial temperature of the Universe has to be high 
enough to assure a sufficient heavy particle production and their kinematic equilibrium. 
It is then quite remarkable that in this case neutrino mixing data favour a simple picture 
of leptogenesis where just decays and inverse decays provide a very good approximation 
[7] and the light neutrino spectrum is hierarchical, with a stringent upper bound on the 
absolute neutrino mass scale rrii < 0.1 eV [8, 9, 10]. The role of the heavy decaying parti- 
cles is played by three new heavy RH neutrinos whose existence is predicted by the seesaw 
mechanism and whose mass matrix is related to the light neutrino mass matrix m v by the 
well known seesaw formula, 

m u = -m D —m T D , (1) 

where mp = v h is the Dirac neutrino mass matrix generated by the Yukawa coupling 
matrix h and v is the Higgs vacuum expectation value. 

It is always possible to work in a basis where the heavy neutrino mass matrix is 
diagonal, M = diag(Mi, M 2 , M 3 ) = D M , with M 1 < M 2 < M 3 . Moreover one can also 
simultaneously diagonalize the light neutrino mass matrix m u by mean of the unitary 
matrix U, such that 

tfm v U* = -D m , (2) 

where D m = diag(mi, m 2 , m 3 ) and mi < m 2 < m 3 . The unitary matrix U can be 
identified with the PMNS mixing matrix in the basis where the charged lepton mass 
matrix is diagonal 1 . In this way the seesaw formula (1) gets specialized in the following 
1 This is rigorously valid in the limit where the heavy neutrinos arc infinitely heavier than the Dirac 
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way: 

D m = U^m D Dl}m T D U\ (3) 
This expression can be also re-casted as an orthogonality condition, 

n n T = n T n = i , (4) 

for the f2 matrix defined as [11] 

Sl = D- 1 '*Vim D D£< 2 (5) 
and whose matrix elements are then simply given by 

0^ = -^= = -^=, (6) 

where rh D = m D and h — Wh. The orthogonal matrix Q is fully determined by three 
complex parameters. Inverting the relation (5) one obtains 

m D = U Dll 2 n D]i 2 . (7) 

Together with the 6 (real) parameters necessary to specify the Q matrix, the seesaw 
parameters are given by the three light neutrino masses, the 3 heavy neutrino masses 
and the 6 parameters (3 mixing angles and 3 phases) that fix the mixing matrix U [12]. 
The expression (7) is particularly enlightening because the 'observable' quantities are 
grouped apart on the RH side, whereas, on the LH side, mo can be regarded as the 
theoretical quantity to be specified within some model or set of assumptions 2 . Among 
the 18 observables, we can actually measure directly only a few of them. A complementary 
useful information can however be inferred when leptogenesis is assumed to explain the 
matter-anti matter asymmetry of the Universe. 

Neutrino mixing data provide two important piece of information on the neutrino mass 
spectrum. In the case of normal (inverted) hierarchy one has 

ml -ml = Am 2 tm (Am s 2 ol ) , (8) 
ml -ml = Am s 2 ol (Am 2 tm ) . (9) 

The third undetermined degree of freedom, the absolute neutrino mass scale, can be 
conveniently expressed in terms of the lightest neutrino mass m\. The two heavier neutrino 

mass matrix eigenvalues and the seesaw formula is exact. 

2 The 3 heavy neutrino masses arc fundamental quantities but, through leptogenesis, they can be also 
regarded as direct 'observable' quantities, in the general meaning explained further on. 
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masses are then given by 

ml = m 2 + m 2 tm , (10) 
m 2 = m 2 + m s 2 ol (Am 2 tm ), (11) 

where we defined m atm = a/ Am atm + Arn^, and m sol = ^Am s 2 ol . The latest measure- 
ments give, combining the results from accelerator and atmospheric neutrino experiments 
[13], 

Am^ tm = (2.3±0.4)xlO- 3 eV 2 , (12) 
and from solar neutrinos experiments [14] 

Am s 2 ol ~ (8.2^| x 1(T 5 ) eV 2 , (13) 

implying 

m atm = (4.9 ± 0.4) x 1(T 2 eV and m sol = (9.0l^ 5 ) x 10" 3 eV . (14) 

Neutrino mixing experiments also measure two of the three mixing angles, 6*12 and #23, and 
place an upper bound on #i 3 . In future they should be able to place, or hopefully measure, 
a stringent constraint on the Dirac phase. On the other hand they are insensitive to the 
two Majorana phases but some information can be extracted from neutrinoless double 
beta decay experiments. 

It could seem that there is no hope to get information on the other 'high energy' pa- 
rameters, the three heavy neutrino masses and the 6 parameters necessary to describe the 
orthogonal Q matrix. However, within leptogenesis, the baryon asymmetry is explained 
as the relic trace of a high energy scale dynamical process, whose final value depends also 
on those seesaw parameters that escape conventional experimental investigation. 

Cosmic Microwave Background observations measure the baryon asymmetry of the 
Universe with a few per cent precision that is further improved when the information 
from large scale structure is included. A recent measurement, expressed in terms of the 
baryon to photon number ratio at the recombination time, is given by the WMAP [15] 
plus SLOAN [16] experiments that find 

jyCMB = 3 ± 0) x 10 -10 _ (15) 

There are different ingredients that enter the calculation of the final baryon asymmetry in 
leptogenesis. In principle the predicted final asymmetry depends on all seesaw parameters, 
however the dependence on many of them is marginal and in this way the information 
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from leptogenesis, that in general would spread out through all parameters, focuses on a 
restricted subset on which useful constraints can be derived. This marginal dependence on 
many of the parameters, as we will see those associated to the asymmetry generation and 
wash-out from the two heavier RH neutrinos, follows not only from intrinsic features of 
leptogenesis but also because neutrino mixing data, the large mixing angles and the values 
of m atm and m so i, select regions in the space of parameters where the final asymmetry 
depends on a limited number of seesaw parameters. Moreover it is remarkable that, 
within leptogenesis, there is a nice conspiracy between the observed baryon asymmetry 
and neutrino mixing data, such that their compatibility is realized for a wide variety of 
models of neutrino masses and mixing, without particular tuning. There is however an 
important limitation represented by a lower limit on the lightest heavy neutrino mass 
Mi [17, 8] and a related one on the initial temperature T[ n of the radiation dominated 
Universe expansion [7] that, within an inflationary picture, can be identified with the 
reheating temperature T rc h [18]. 

In Section 2 we outline the different steps and approximations that lead to predictions 
of the baryon asymmetry depending only on a restricted number of seesaw parameters. 
This will enlighten the importance of the maximum CP asymmetry and of the effective 
leptogenesis phase to describe such predictions. We will concentrate our attention on 
these two quantities in the subsequent analysis. A useful tool is provided by the 'seesaw 
geometry' that we introduce in Section 3 and consists in representing the orthogonal Q 
matrix in the complex plane. In this way, in Section 4, we are able to determine the 
leptogenesis CP bound, generalizing and clarifying existing results [17, 10, 19] and speci- 
fying the conditions on the Q matrix for the effective leptogenesis phase to be maximal. 
In Section 5 we find an expression for the effective leptogenesis phase placing an interest- 
ing lower bound and studying the conditions on the Q matrix parameters for successful 
leptogenesis. Our procedure is phenomenological, without any assumption on the matrix 
of the Yukawa couplings, while the compatibility between different theoretical predictions 
onmj) with the observed baryon asymmetry and neutrino mixing data has been investi- 
gated in many works [20] . On the other hand, we will often point out different interesting 
connections between our results and neutrino models, in particular in connection with 
the lower bounds on M 1 and on T Teh . In Section 6 we discuss a new scenario of thermal 
leptogenesis where the final asymmetry is generated by the decays of the second lightest 
RH neutrino. This is possible when the wash-out from the lightest RH neutrinos is weak 
and we will see how this new possibility solves the well known problems of this regime, 
namely the dependence on the initial conditions and the fine-tuning of the effective neu- 
trino mass rh\. In Section 7 we examine the conditions of validity of the assumption of 
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a hierarchical heavy neutrino spectrum. In Section 8 we summarize and draw some final 
conclusions. 



2 Seesaw parameters and leptogenesis 

The final baryon asymmetry from leptogenesis is, in general, a function of all seesaw 
parameters. More explicitly one has 

Vb = ?7s(mi, Mi, 0i3, 5, if!, if 2 , ^; m atm , m soh 9 12 , 2 3) , (16) 

where the parameters measured in neutrino mixing experiments are indicated after the 
semi-colon. Their values justify some approximations and simplifications, such that the 
final asymmetry depends on a restricted number of parameters [8]. 

An important approximation is to neglect the asymmetry generated by the decays of 
the two heavier RH neutrinos calculating the final baryon asymmetry only from the decays 
of the lightest ones. We will consider in Section 6 a special scenario where this approxi- 
mation does not hold while in Section 7 we will specify the conditions for its validity. On 
general grounds this approximation is well justified because neutrino mixing data favor a 
picture where, even for a mild heavy neutrino mass hierarchy, the asymmetry generated 
by the two heavier neutrinos is efficiently washed-out and can be safely neglected. This is 
true for two reasons. The first is that large mixing angles produce a flavor misalignment 
such that the asymmetry produced from the decays of the two heavier RH neutrinos can 
be washed-out mainly by the inverse decays of lightest RH neutrinos and by other sec- 
ondary processes. The second is that, the experimental fact that m atm ,m so i 10~ 3 eV, 
makes typically such a wash-out highly efficient and only what is produced by the lightest 
RH neutrino decays around the baryogenesis temperature Tb <C M\ survives until the 
present. This also implies that one has to assume the initial temperature to be slightly 
higher than Tg [8, 7]. With this assumption, one can more generally neglect any ini- 
tial asymmetry and assume an initial thermal abundance of the lightest RH neutrinos 
without worrying about the exact initial conditions and of a detailed description of the 
thermalization process [10, 7]. In this case tjb from leptogenesis can be calculated as 

7 

The sphaleron converting coefficient a sp h ~ 1/3 is the fraction of the B — L asymmetry 
that ends up as a baryon asymmetry because of sphaleron conversion [21]. The number 
of photons at recombination (in the portion of comoving volume containing, on average, 
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one Ni in ultra-relativistic equilibrium) can be calculated assuming a standard thermal 
history and in this case i\^ ec = 4 (g SM + 7/4)/(3 # rcc ) ~ 37. 
The CP asymmetry parameter is defined 3 as 

r - f 

£l = -fTT- (18) 

where T (T) is the decay rate into leptons (anti-leptons) and it gives the (B—L) asymmetry 
generated per single decay on average. 

The final efficiency factor K{ takes into account both the heavy neutrino production 
and the wash-out due to different processes, in particular inverse decays. It is defined in a 
way to be one in the asymptotic limit of vanishing decaying rate, initial ultra-relativistic 
thermal neutrino abundance. If the condition Mi ^\ mf <C 10 14 GeVm 2 tm holds, then 
the efficiency factor depends only on the effective neutrino mass mi, defined as [22] 

^ (4^)n (19) 
Mi v ' 

Plugging the relation (7) into the effective neutrino mass definition, one easily gets [23] 

ffii = mipn +m 2 p2i +m 3 p 31 , (20) 
where we have introduced the polar representation 

PU < '"■■■ (21) 

with pij = > 0. Since the first row of the Q matrix will play a special role in our 
discussion, it is useful to simplify the notation such that pji — > pj and ipji — > ipj. The 
orthogonality of the Q matrix implies that 

ft* + fi 2 2 i + ^ = 1 , (22) 

from which it follows that rh\ > mi [23]. This is the only fully model independent 
restriction on rh\. For configurations such that 

E Pi ~ E n n = 1 ( 23 ) 

3 3 

one has rh\ < m 3 . Models with rh\ ^> m 3 rely on the possibility of fine tuned phase 
cancellations, since only in this way one can have p 3 ^> 1, with p 2 or pi, or both, com- 
parable to p 3 . Independently on whether these models can provide or not a satisfactory 
3 To have it positive, we have introduced a minus sign compared to the usual definition. 
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solution to describe neutrino masses and mixing and whether they can or cannot be jus- 
tified within some general theory, one has also to consider the leptogenesis upper bound 
mi < (0.009 eV) (Mi/10 10 GeV) - 8 , valid for M 1 < 10 14 GeV 4 . 

Conversely, there are models with pi ~ 1 ^> Pj^i- Since the values of the Re(f2^)'s 
determine the contribution of the heavy neutrino Nj to the determination of the light 
neutrino mass rrii [12, 24], these kind of models are characterized by the dominance of 
the lightest RH neutrino in the determination of the lightest neutrino mass In the 
strict limit Pi = 1 and pj^ = 0, one has necessarily rhi = rrii < m 3 . If the other two light 
neutrino masses, rrij^i, are also dominantly determined, by one of the two heavier RH 
neutrinos, and if Im(fi?-) vanish, then one has a particularly simple sub-set of the class 
models with sequential RH neutrino dominance [25]. In the limit case of exact dominance, 
there are 6 possible choices for the Q matrix given by, 
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plus the three cases obtained by exchanging the second and third row. Obviously for these 
exact limits, the CP asymmetry and the effective leptogenesis phase vanish. However, it 
is enough to add small complex perturbations to have non vanishing CP asymmetry and 
we will see that this can be even maximal if one starts from the first matrix in Eq. (24). 
An important issue for leptogenesis is whether rh\ is larger or smaller than the equilibrium 
neutrino mass m* ~ 10~ 3 eV, that sets the transition between the weak (for rhi < m*) 
and the strong wash-out regime (for rhi ^> m*). In the latter case thermal leptogenesis 
can work at its best with no dependence on the initial conditions and minimal theoretical 
uncertainties [7]. The possibility to have m x < m* relies on the fulfilment of a few 
conditions. First of all, since rhi > m x , one must have mi < 10 _3 eV m so i <^ m atm 
and thus it requires fully hierarchical neutrinos. From Eq. (20) it follows that p2,P3 ^ 
m*/m so i, m^/m atm <^ 1 and then, because of the orthogonality of Q, that pi ~ 1. This 

4 This bound can be simply obtained inverting the lower bound on Mi given in Section 5 (cf. (61) 
and (65)). For M\ > 10 14 GeV the bound is even more restrictive. The maximum possible value for mi, 
independently on the Mi value, is obtained at Mi ~ 5 x 10 13 GeV and is about rhi ~ 10 eV, anyway 
much higher than the naive expectation, often quoted in the literature, mi < m 4 ~ 10 -3 eV, from the 
full out-of-equilibrium condition. This bound holds under the conditions on e x studied in Section 7. 
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selects models with Q matrices that are perturbations of the complex 23-rotation 



ft = 



( 1 
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and that can be obtained performing additional 'small' rotations in the planes 13 and 
12. A specific realization has been shown in [10]. For these models one has m x ~ 
v 2 \h\i\/Mi and thus the requirement mi < rhi < m* implies \hn\ < a = Mim^/v 2 ~ 
10~ 3 a/ Mi/10 10 GeV. At the same time one has to require p2(3) <C ^*/ m soi(atm) — 
0.125(0.02), for normal hierarchy, and this is equivalent to have \hn\ <C a. On the 
other hand in the limit rhi — > mi, as we will see, S\ vanishes like 1 — mi /% and thus one 
has to impose also rhi ^> mi in order to maximize the asymmetry. In this way one arrives 
at the conditions \hu\ <C \h 2 i\ + \hsi\ < a. There are no compelling reasons to exclude 
these kinds of models, but it is clear that they require a good amount of fine tuning to be 
realized. Note also that the request that the lightest neutrino mass is dominantly deter- 
mined by the lightest RH neutrino mass is such that these models are excluded when one 
imposes the experimental requirement of large mixing angles and the theoretical require- 
ment that the neutrino Yukawa matrix resembles hierarchical quark matrices [25] . These 
difficulties are in addition to the problem, within leptogenesis, of the dependence of the 
final asymmetry on the initial conditions. However, in Section 6 we will discuss a new 
scenario where rhi % but the problems of the fine-tuning of rhi and of the dependence 
on the initial conditions are solved. For the time being we will refer to the more general 
case of strong wash-out regime for rhi m*- 

Let us now consider the CP asymmetry. A perturbative calculation from the interfer- 
ence between tree level and vertex plus self energy one-loop diagrams yields [26] 



£i = 



Im [(fct h )l] 



8tt ^ 

i=2,3 



X 



f | M H + f ( M * 



The function fy, describing the vertex contribution, is given by 

1 + x 



M 2 



(26) 



fv(x) 



1 - (1 + x) In 



x 



while the function f s , describing the self-energy contribution, is given by 

y/x 



fs(x) 



1-x 



(27) 



(28) 



Notice that this last expression is valid only when the two mass differences Mj — M 1 
are larger than the difference of the respective two decay rates. In the limit x 3> 1, 



9 



corresponding to have a RH neutrinos mass hierarchy with Mf S> Mf , one has 

/ V (x)+/ fl (x)~-JL. (29) 

In this limit the expression (26) simplifies into [27] 

3 M l Im [(frUM-^/Qn] 

£l " T&7" (/^ (3U) 

and in Section 7 we will study the conditions of its validity. Replacing h with Q (cf. (6)), 
one then obtains [17] 



ei ~ ^(Mi.mi.mi.Q !) = — = /3(mi, rai, Q J , (31) 

J l07T f z J 

where we have introduced the convenient dimensionless quantity 

*™w=ibt£m- (32) 

We have now considered all the different ingredients that enter the prediction of the final 
baryon asymmetry within the approximation of a hierarchical heavy neutrino spectrum 
and from Eq. (17) we can see that r\ B = ^(M^ mi, mi, f^) . 

Notice that if the value of rh\ is fixed, then a generic configuration of three because 
of the orthogonality condition (22), depends only on 3 independent parameters. Therefore, 
the final asymmetry depends just on 6 unknown seesaw parameters. This is true within 
the validity of the approximation of a hierarchical neutrino spectrum that made possible 
to cancel out a dependence on M 2 , M 3 and on two parameters of the Q matrix. Except 
for the absolute neutrino mass scale m 1; the other 5 parameters cannot be measured 
in terrestrial experiments. In Section 6 we will study a situation where, even though 
the heavy neutrino spectrum is still assumed to be hierarchical, the dependence of the 
final asymmetry on the second lightest RH neutrino cannot be neglected and additional 
parameters enter the prediction of the final asymmetry. Conversely, in Section 7, we will 
study the conditions for the validity of the assumption of hierarchical heavy neutrino 
spectrum and where this breaks down, such that, again, additional parameters have to 
be taken into account for the determination of the final asymmetry. 

The efficiency factor Kf depends only on fhi, mi and M 1 , implying that the dependence 
of the final asymmetry on the f^i's arises only from E\. It is then quite interesting to 
search for those particular configurations of the Jl^'s such that, for fixed values of m\ 
and mi, the CP asymmetry is maximum and, for this purpose, it is useful to define an 
effective leptogenesis phase 5l [28] such that 

/3(mi,mi,fiy = Anax(mi,mi) sin5 L (mi, fh u Q 2 ^) . (33) 
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The function /3 max (mi,mi) > is the maximum value of (5 corresponding, when plugged 
into Eq. (31), to a maximum CP asymmetry e^ ax (M 1 ,m 1 ,fhi) > 0, realized for those 
particular configurations with sin^ = 1. Notice that a trivial bound on the CP asymme- 
try is, by definition, E\ < 1. It is however possible to find a more restrictive non trivial 
bound within the validity of Eq. (30) [28, 17, 10, 19]. The maximum baryon asymmetry 
is given by 



where d = a sp h/A^ cc and it depends only on the three seesaw parameters mi, m x and 
Mi. Successful leptogenesis implies ?7g iax (M 1 , m 1? rhi) > r]g MB and this yields interesting 
leptogenesis constraints on mi, mi and M x [17, 8] and in particular an upper bound on 
the absolute neutrino mass scale mi [8, 9, 10] testable in cosmological, neutrinoless double 
beta decay and Tritium beta decay experiments. 

For an accurate evaluation of these constraints, it is thus necessary to determine the 
function /5 max (mi,mi) and it is important to understand all properties to be fulfilled in 
order to saturate the bound. This means a determination of the class of configurations of 
f2j\'s such that sin 5l = 1- At the same time we will be interested in determining the phase 
sin#L for an arbitrary configuration of 1X,i's and its connections with neutrino models. In 
this way one can answer interesting questions on the stability of the CP bound, that means 
whether small variations of the seesaw parameters determine a small or a large variation 
of the leptogenesis phase. This will make possible to realize whether the condition of 
maximal phase should be regarded as a reasonable or a fine tuned assumption. Moreover, 
placing a lower bound on the same leptogenesis phase will provide interesting information 
on the geometrical parameters. 

3 Seesaw geometry 

Let us represent the three il^'s in the complex plane. The orthogonality condition fixes 
the sum of the three to start from the origin and to end up onto the real axis at the point 



j 

as shown in Fig. 1 for a generic configuration (dashed line arrows). In addition to the 
polar representation (cf. (21)), it is also useful to indicate explicitly the real and the 



^ iax (Mi,mi,mi) = d£5 nax (Mi,mi,m 1 )K f (Mi,m 1 ,m 1 ), 



(34) 
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Figure 1: seesaw geometry. An arbitrary Ctj-y configuration with sign(Ya) = — sign(Y) (dashed 
line arrows) and its 'dual' one with sign(>3) = sign(Y) > (dotted line arrows). The region 
inside the short-dashed line is that one spanned by the areas delimited by the configurations 
with = Pj < 1 and the segment [0,1]. Configurations corresponding to models with 

Pi ~ 1 > pj^i lie within a region, schematically indicated with the box, around the [0, 1] 
segment on the real axis. 



If one defines that class of models such that pj < 1, then the area delimited by the three 
fi^'s and the segment [0, 1] is confined to lie within a well denned region in the complex 
plane inside the dashed line in Fig. 1. For these models one has (cf. (20)) 

~ I ?7iatm (fully hierarchical neutrinos) 

mi<mi+m 2 + m 3 ~ (37) 

I 6 mi (quasi-degenerate neutrinos) 

Models with pi ~ 1 3> Pj^i correspond to all those geometrical configurations with small 
deviations from the segment [0, 1] along the real axis (see the box in Fig. 1). In particular, 
as we have already seen in the previous section, only for models with p\ ~ 1 p<i-,Pz (cf. 
(25)) one can have rh\ < m*. 

Conversely, the possibility of having rh\ ^> m 3 relies on models with strong phase 
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cancellations with no dominance of just one pj, while at least two of them have to be 
comparable and much larger than one. The two (dual) configurations depicted in Fig. 1, 
are an example of such a possibility. 

An interesting aspect of this geometrical representation is a comparison with the CP 
violation in the weak interactions of quarks, described by the CKM matrix, and in neutrino 
mixing, described by the PMNS matrix. In the quark sector a geometrical representation 
of CP violation and of the CKM matrix is provided by the unitarity triangles, whose area 
is given by the Jarlskog parameter [29]. If the triangles squeeze to a segment, then CP 
violation vanishes. This is an easy way to understand why a non vanishing CP violation in 
the quark sector and in neutrino mixing requires at least three generations. On the other 
hand in the case of the CP violation governed by the orthogonal Q matrix two generations 
are enough, since instead of a triangle one has, in general, a quadrilateral where one side, 
the real segment [0, 1], is present anyway. Thus, within the seesaw mechanism, one has a 
kind of CP violation that is different from the usual one in the quark sector (or from the 
possible one in neutrino mixing), since in the first case this is described by an orthogonal 
matrix, while in the second case it is described by a unitary matrix. 



4 CP asymmetry bound 

Let us now study the leptogenesis CP asymmetry bound taking advantage of the 'seesaw 
geometry'. Using the rh\ definition (cf. (20)) and the orthogonality condition (cf. (22)), 
the function (3(m,i, rhi, ^ 2 J, defined by the Eq. (32), can be re-casted as 

P( mi , m,, ft 2 = 32 = 21 , 38 

J m atm m 1 

with Y = Y 2 + Y 3 . Notice that from any configuration such that sign(F) = — sign(Y 3 ), 
one can always obtain two configurations with equal signs, sign(F) = sign(Y 3 ) = ±1, 
same value of rh\ and such that \(3(mi, mi, ^ 2 i)| is higher. Since the observed baryon 
asymmetry is conventionally defined to be positive, then the CP asymmetry, the function 
j3 and the effective leptogenesis phase 8 L have to be positive too and thus we will always 
refer, in the following, to configurations with positive Y and Y 3 (an example is shown in 
Fig. 1). 

It is easy to find the upper bound /3 max (mi, mi) (cf. (33)) with mi free to change while 
mi is fixed. If we write rhi as 

mi = mi y/Y* + (1 - Xf + m 2 y/x[+Y? + m 3 y/x$ + Y? , (39) 
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then (3(mi, mi, Q^) is maximized in the limits Y 3 ^> X 3 , Y 2 ^> X 2 and Y ^> 1, implying 
mi ^> m 1? and it gets further maximized for Y 2 = 0. In this way one obtains the upper 
bound [17] l3 mSLX (mi , mi) = (m 3 - mi)/m atm , corresponding to (cf. (31)) 

m 3 — mi 



c max 

(Mi) 



Wlatm 



(40) 



where we defined 



3 M^^ 10 _ 

l07T v z 



Mi 



10 10 GeV 



f-^ 2 -) , (41) 
V0.05eV/ ' V ; 



giving the absolute maximum of the CP asymmetry in the limit of hierarchical neutrinos 
[28, 17] for mi/m atm — > 0. For a given finite value of rh\ there is a more restrictive upper 
bound that in general, introducing a function /(mi, mi) < 1, can be written as [10] 



/3max(rai,rai) 



m 3 — mi 

^atm 



/(mi, mi) < 1, 



corresponding to 



e^Mi^i.mi,^) < £ max (Mi 



m 3 — mi 



/(mi, mi) . 



(42) 



(43) 



^atm 

For rhi/mi — > oo one has /(mi, oo) = 1 and the limit (40) is recovered. 

We have now to calculate the function /(mi, mi) for finite values of rh\. Before dealing 
with the general case, for arbitrary mi, let us consider the two interesting asymptotic lim- 
its of fully hierarchical neutrinos, for (mi/m so i) 2 — > 0, and of quasi- degenerate neutrinos, 
for (mi/m atm ) 2 -> oo. 



4.1 Fully hierarchical neutrinos 

In this case /3 max = 1 and there is no global suppression. Moreover one has 

fhi =m 2 p2 + m 3 p 3 , 



(44) 



with m 3 = m atm and m 2 = m so i {\fm 2 &tm — m 2 ol ) in the case of normal (inverted) hierarchy. 
This implies that, for any change of configuration such that p 2 and p 3 remain constant, 
the quantity rhi is constant too, while p x can be arbitrarily modified. Hence, for a given 
value of rh\, one has that /9(0, rh\, Q 2 X ) is maximum for a configuration with X 2 = X 3 = 0, 
such that the Eq. (44) gets specialized into rh\ = m 2 Y 2 + m 3 Y 3 . Using this expression 
one can replace Y 3 into the Eq. (38), getting 



Anax(0,mi) = /(0,mi) = max 



m 2 I m 2 



mi 



m 3 



(45) 
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that is clearly maximized for Y 2 = 0, corresponding to Y 3 = Y and rh\ = m 3 Y 3 . Therefore, 
one has very simply f(m 1 = 0, mi) = 1 showing that the case m 1 = yields, for a fixed 
Mi, an absolute maximum of the CP asymmetry given by e max (Mi) (cf. (41)), as we 
have already found as a particular case of mi/mi > 1. In Fig. 2 we show a generic 
configuration with a given fhi (solid line arrows) and the corresponding configuration 
X3 = Y3 = Y2 = 0, with the same value of mi, for which the CP asymmetry is maximum 
(dashed line arrows). It is worth to notice that this class of configurations correspond to 
the case when the lightest neutrino mass is dominated by the lightest RH neutrino mass 
and thus for any other kind of model there will be necessarily some phase suppression 5 . 

4.2 Quasi-degenerate neutrinos 

In the quasi degenerate limit, for (mi/m atm ) 2 3> 1, one has mi ~ m 2 — m 3 and the 
expression for rh x (cf. (20)) becomes simply rhi ~ mi £, with £ = ^\ pj. Therefore, 
the condition m x = const is equivalent to select all those configurations for which £ is 
constant. Then it simply turns out that /3(mi, mi, Jl^) is maximum for a configuration 
such that Y 2 — X 2 — and X 3 = 1/2, as shown in Fig. 2 (dashed line arrows), for the 
simple reason that only for this configuration both Y 3 and Y are simultaneously equal to 
their maximum value. This is easy to be understood in a geometrical way: it corresponds 
to 'stretch' the string with length £ as high as possible in order to maximize Y, while 
at the same time the value of Y3 is maximized by having X 2 — Y 2 — 0. Since in the 
quasi-degenerate limit (neglecting terms C , (m atm /'^i) 2 ), one has m^ tm ~2mi (m 3 — mi), 
it is easy to obtain the result (cf. (38) and (42)) 



already obtained in [19] with a different procedure and using the approximation m so \ = 0. 

5 In [30] it was found that, within single RH neutrino dominated models, where the largest neutrino 
mass TO3 is dominated by one RH neutrino, the CP asymmetry e\ is larger when the largest neutrino mass 
TO3 is dominated by the heaviest RH neutrino compared to when is dominated by the lightest. This can 
be regarded as a particular case of our more general result: if the lightest RH neutrino dominates to 3 then 
it cannot dominate mi , because of the f2 orthogonality, and thus there must be a phase suppression. On 
the other hand having m 3 dominated by the heaviest RH neutrino is compatible with mi dominated by 
the lightest RH neutrino. The point is that the RH neutrino dominance directly relevant for leptogenesis 
is that of mi and not that of m 3 . 




(46) 
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X=Re(£$ 1 

Figure 2: An arbitrary configuration for a given rh\ value (solid line arrows) and the associated 
ones that, for the same rh\ value, maximize the CP asymmetry in the limits of fully hierarchical 
neutrinos (dashed line arrows) and degenerate neutrinos (dot-dashed line arrows). 

4.3 The general case 

Let us now find the configuration that maximizes the CP asymmetry and calculate 
/(mi, mi) for arbitrary values of m\. Notice that in both the two asymptotic limits 
we have found VL% X = and so one can wonder whether this condition holds in general for 
any value of m\. This can be proved showing that starting from a generic configuration 
with p 2 7^ one can always find a configuration with the same value of mi, such that p 2 
is lower and (3 is higher. It is useful to split the demonstration in two steps. 

The first step consists in showing that starting from a generic configuration one can 
always find a configuration, with higher (5 and same value of rh\ and such that Y 2 = 
(i.e. Y — y 3 ). The procedure to build such a configuration is depicted in Fig. 3 and in 
Fig. 4, while we address to the Appendix A for more details. The second step consists in 
showing that within the class of configurations with a given value of rh\ and with Y 2 = 0, 
the value of (3 is further maximized when X 2 = and thus, in conclusion, when f^i = 0. 
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Y 




Figure 3: Procedure that brings from a generic configuration with Y 2 > (dashed arrows) or 
Y 2 < (dotted arrows) to one with Y 2 = (dash-dotted arrows) and from this to one with 
= (solid arrows) such that mi does not change and the CP asymmetry increases. 



Looking at the two configurations in Fig. 3 and in Fig. 4 with Y 2 = (dash-dotted line 
arrows), it is quite intuitive geometrically that Y increases when X 2 decreases and this is 
confirmed by an analytic calculation. Indeed, the expression (38) simplifies into 

P(m 1 ,m 1 ,tf jl ) = ^Y, (47) 
J mi 

where now Y = Y(rhi, X3, X 2 ) is a function of just three parameters. Since rh\ is fixed, 
then (3 is maximum when Y is maximum. In the quasi- degenerate case it is easy to 
understand that Y is maximum when X3 = (1 — X 2 )/2 and is further maximized when 
X 2 = such that X 3 = 1/2, confirming the result that we previously obtained (cf. 
subsection 3.2) in another way. In the hierarchical case (mi = 0) it is simple to obtain 
from the Eq. (44) that 

m 3 m 3 

that is maximum for X 3 = X 2 = 0, confirming again what we have previously found 
in subsection 3.1. In the general case, starting from the definition of fh\ (cf. (39)) 
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Figure 4: An example of the procedure that brings to a configuration with maximum asymmetry 
starting from one with Y2 > (dashed arrows). In this example by decreasing P2 the config- 
uration gets first transformed into one with 922 = ^3 (dotted arrows). From this recursively, 
as explained in the Appendix, one still ends up with a configuration with Y2 = (dot-dashed 
arrows) and from this again to one with f^i = (solid arrows). 

with Y = Y 3 and taking the derivative with respect to P2, one can easily show that 
dY/dp 2 < 0. Thus the maximum value of Y is obtained when Q^i = 0. Two examples of 
such configurations are shown in Fig. 3 and in Fig. 4 (solid arrows). We can thus conclude 
that for a given value of mi and rh\, the CP asymmetry is maximized for configurations 
with f^i — 0- 

For a given value of fhi, these configurations are described just by one parameter, and 
we can conveniently choose X — 1 — X\ — X3. Therefore, the expression (39) for rh\ can 
be re-casted as 

VV 2 + (1 - X) 2 + m 3 VY 2 + X 2 . (49) 

The last step is thus to maximize the CP asymmetry and this means to find the maximum 
value of Y entering the expression for (3 (cf. (47)), with respect to X. In this way one can 
obtain the CP asymmetry bound for given, arbitrary, values of fhi and mi, not depending 
on the geometrical parameters Vt 2 v This can be found just taking the derivative of the 
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expression (49) with respect to X and imposing (dY/dX) x =x i , = 0. Doing this, one 
arrives to the simple relation m 3 cosy? 3 = m-i cos</?i, or explicitly in terms of X+ and Y m 

m 3 X^ m 1 (l-X ir ) 

(50) 



This expression, together with the Eq. (49) for fhi, defines the functions X±{mi,rhi) and 
7 m (mi,mi), such that (3 and thus the CP asymmetry are maximum. Notice that the 
function f(mi,fhi) (cf. Eq. (42)) is related to F m (mi,mi) by 

/(mi, mi) = m3 + mi y m ( mii mi) . (51) 
mi 

In Fig. 5 we have plotted / and X* vs. x = mi/ fhi and for different values of mi/m atm - 
It is also possible to derive explicitly analytic expressions for the functions X* and Y m . 
However these are algebraically quite involved and not useful to be written down, while 
the two equations (49) and (50) provide a much more useful physical insight and at the 
same time simple and useful explicit analytical expressions can be easily derived in some 
interesting limits: 

• in the limit of fully hierarchical neutrinos one obtains X* = and Y m = mi/m 3 , 
implying the already well known result /(mi, mi) = 1; 

• in the limit of quasi- degenerate neutrinos, for mi ~ m 3 , one gets from the Eq. 
(50) X* = 1/2, that inserted in the Eq. (49) gives once more /(mi, fhi) = 
a/1 - (mi/mi) 2 ; 

• in general, since < mi/m 3 < 1, one has < X* < 1/2; 

• if fhi — mi then X = Y = f = and vice versa; 

• in the limit fhi ^> mi one gets easily Y m (mi,fhi ^> mi) = fh\j (mi + m 3 ), corre- 
sponding to the already known result / '(mi, fh\ mi) = 1, and also X*(mi,mi ^> 
mi) = mij '(mi + m 3 ). 

It is also easy to go beyond the limit of fully hierarchical neutrinos and find an expression 
for the function / '(mi, fhi) valid in the regime mi/m atm < 1, corresponding to m x < 
m 3 /\/2- The equation (50) can be conveniently re-casted in the following form 

XjY m mi (1-X*)/Y m 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x=m 1 /m 1 

Figure 5: function / calculated from the Eq.'s (49), (50) and (51) and plotted vs. x = mi/fhi 
for TOi/m a tm = 5, 2, 1, 0.5, 0.2, 0.1 (from the top to the bottom solid line). The hierarchical limit 
for / (cf. (56)) is plotted for mi/m a t m = 1, 0.5, 0.2, 0.1, 0.01 (from the top to the bottom dashed 
line), while the quasi-degenerate limit (cf. Eq. (46)) corresponds to the dot-dashed line. 



from which one can see that 



X* mi 



Ym y'mj — m\ 



(53) 



and from this it follows that X+ < mi / \V2 m 3 ) . Neglecting terms 0{m\/m\) one can 
then find the approximate relation 

mi Y m 



m 3 v/TTTT 

while the Eq. (49) becomes simply 

mi a/F^ + 1 + m 3 Y m = rhi ■ 
It is in the end easy to obtain the expression 



/(mi, mi) 



m s — mi \ 1 + 



m3 — mi 



(54) 



(55) 



(56) 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x=m 1 /m 1 



Figure 6: X+ calculated from the Eq.'s (49) and (50) as a function of x = mi/mi and for the 
indicated values of mi/m atm . 

valid at the first order in m 1 /m 3 , while at the zero-th order in mi/m 3 and mi/m 3 (thus 
for m atm ~ rri3 > mi > mi) it becomes simply /(mi, mi) = 1 — m\jm\. The expression 
(56) has been found in [10] assuming X* = and used to get an upper bound on the 
absolute neutrino mass scale m^ < 0.1 eV. This function clearly does not describe the 
function /(mi, mi) for any value of mi. In Fig. 5 one can see (dashed lines) how for 
increasing values of mi/m atm , from 0.01 to 1, the function saturates to its asymptotic 
limit / = 1 — (mi/fhi) 2 and so, in particular, it does not reproduce the correct bound 
represented by the Eq. (46) [19]. However for values m-i ~ 2m atm ~ 0.1 eV, those 
corresponding to the bound on neutrino masses, it underestimates the bound on the CP 
asymmetry at most of ~ 20% (see Fig. 5 and compare with the solid line for mi/m atm = 2) 
and the neutrino mass upper bound itself is therefore underestimated only of ~ 5% 
(corresponding to ~ 0.01 eV) [7], much below the precision of the bound itself. This, 
together with the result that we have just shown, for which this function is correct in the 
limit (mi/m^) 2 < 1, fully justifies its approximate use for an arbitrary value of mi and 
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is particularly interesting in connection with a search of the neutrino mass upper bound 
when some restriction on the space of relevant parameters is introduced, such that the 
neutrino mass bound becomes more stringent and thus it falls in the region (mi/m 3 ) 2 < 1. 
This is for example the case when a cut-off < 10 12 GeV on the value of M\ is imposed 
[10, 32]. In Section 6 we will also show another ueful application of the Eq. (56). We 
have also compared the exact /(mi, mi), given by the equations (49), (50) and (51) and 
plotted in Fig. 5 (solid lines), with the results obtained in [19] using the approximation 
m &o\ — 0. We have found a slight difference, within a few per cent, for mi ~ fhi, getting 
completely negligible as soon as rh\ is appreciably larger than mi. 



5 The effective leptogenesis phase 

In the previous section we have characterized which are the configurations that maximize 
the CP asymmetry, or equivalently for which the effective phase sinJ^ = 1. Here we want 
to study how easily this condition can be realized and which are the typical values of the 
phase for different classes of configurations. In general we can write, from its definition 
(cf. (33),(38),(42)), 

sin^m^mi,^) = = mi f . + ™l (Y 3 + a 2 Y 2 ), (57) 
J mi /(mi, mi) 

where a = ^/Am^/m atm is given by m sol /m atm or by a/1 - m s 2 ol /m 2 tm , depending on 
whether the case of normal or inverted hierarchy is considered. The expression (57) can 
be also recasted as 

. xf ~ (~\2 \ ™i+ m 3 ( Ps sin y? 3 + <J 2 p 2 sin ^ 2 \ 

sm^m^mi,^) = — — ■ ■ , (58) 

J /(mi, mi) \m 1 p 1 +m 2 p 2 + m 3 p 3 J 

showing that the effective phase reduces to genuine phases in particular cases, while in 
general it is given by a linear combination of them, analogously to how the effective 
neutrino mass rh\ is a linear combination of the neutrino masses, with coefficients given 
by the p/s in both cases. It should be stressed that, in the general framework we are 
considering, it does not depend on the phases that could be responsible for CP violation 
at low energies and contained in the mixing matrix U. The requirement of successful 
leptogenesis implies the lower bound 

^CMB 

sin^m^mi,!! 2 !) > J* — , (59) 

j 77™^ (Mi, mi, mi) 

that, together with the lower bounds on Mi and on T re h, represents another useful way 
to get information on the f2-parameters through leptogenesis. For practical purposes it is 
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interesting to study the effective phase in the two asymptotic limits of fully hierarchical 
neutrinos and quasi- degenerate neutrinos. 



5.1 Fully hierarchical neutrinos 

In this case one has m x = 0, m 3 = m atm and / = 1 6 and the general expression (58) 
becomes 

sin 5 L = \ps sin ip 3 + cr 2 p 2 sin ip 2 ] , (60) 

mi L J 

where now simply rh\ = m 2 p2 + ^atm Pz an d m 2 = m so i (~ m atm ) in the case of normal 

(inverted) hierarchy. The value of the phase in the full hierarchical case is relevant in 

connection with the (3 a) lower bound on Mi given by 

M 1 > 4 - 2 *f G f (M 1 «10»GeV), (61) 

that generalizes the case of maximal phase considered in [8]. When thermal effects on 
the Higgs mass are taken into account and a proper subtraction procedure of on shell 
AL = 2 processes is employed [33], the efficiency factor can be approximately calculated 
accounting only for decays and inverse decays [7]. In any case this provides a conservative 
upper bound, since other effects, in particular those arising from spectator processes [34], 
can only lower the final asymmetry. Then a very good analytical approximation in the 
case of an initial thermal abundance, that gives quite a conservative upper limit in the 
weak wash-out regime, is 

where the decay parameter K is related to the effective neutrino mass through the relation 

K— — ~ 926 — , (63) 
m+ ev 

while the baryogenesis z = M 1 /T value is calculated using 7 

z B (K) -9/8 + 2.35 [\n(l + K)} - 8 . (64) 

A comparison between the analytical approximation (cf. (62) and (64)) and the numerical 
solution is shown in Fig. 7 while the lower bound on M\ is plotted in Fig. 8. 

6 Notice that then simply one has sinJ^ = (3 since /3 max = 1- 

7 This is a different expression from the two ones, practically equivalent, given in [7] and in [32] and is 
such to cure the error of the expression (62) around values K ~ 1, yielding an analytical expression that 
well agrees with the numerical calculation within 10% for all values of rh\. 
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Figure 7: Efficiency factor: the analytic approximation (dashed line) given by the eq.'s (62) 
and (64) is compared with the numerical solution (solid line) of the kinetic equations when the 
wash-out from scatterings is neglected and an initial thermal N\ abundance is assumed [7] . The 
dot-dashed line is the power-law fit eq. (65). 

For 0.5 m so i < rh\ < 2 m atm the expression (62) can be approximated with the power 
law [32] 

0.5 „ _ __ 2 /10~ 2 eV x L2 



fcf(mi) ~ -r ~ 3.5 x 10~ 2 — (65) 

K \ mi J 

and in the range rh\ = [m so i, m atm ] one gets, from the Eq. (61), Mi > (10 10 — 10 11 ) GeV. 

Notice that only for K > 5 the lower bound on Mi does not depend on the initial 

conditions and in particular on a specific description of Aq's production. The lower bound 

at K ~ 5, Mi > 5 x 10 9 GeV, is thus particularly meaningful, since it can be regarded 

as the lowest allowed Mi value without the need of additional assumptions beyond the 

given ones. 

Within the minimal supersymmetric standard model (MSSM) the CP asymmetry and 
its maximum value get double and this implies a relaxation of the lower bound Eq. (61) 
at low rhi, such that [32] 

2.25X10* GeV <<; lQli 

Kf{mi) sxhol 

The expression (62) for the efficiency factor is still valid but the stronger wash-out pro- 
cesses make the equilibrium neutrino mass a bit lower than in the SM case, rn+ ISSM ~ 
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Figure 8: Lower bounds (3a) on Mi for different value of the effective phase in the case of normal 
hierarchy. The solid line corresponds to maximal phase in the case p2 = (cf. (61)). The dashed 
line corresponds to the case /O3 = and is given by the Eq. (61) with svd.5l = a ~ 1/6. The dot- 
dashed line is the intermediate case P3 = a 2 P2 with an effective phase sin^ = 2 a ~ 1/3. The 
box indicates a transition range of rh\ values in the scenario where the asymmetry is generated 
by the iVys decays for p% = pi = 0: for rh\ < m* (K < 1) there is no lower bound on Mi, while 
for fhi > 5 m* (K > 5) there is no allowed value. 



0.0008 eV, such that K ~ 1270mi/eV. In this way the efficiency factor in the strong 
wash out regime gets also lower and the power law (65) becomes [32] 

, ( ^.,4 X1 0-(12JX)'- 2 , (67) 

while the 3<r lower bound on Mi in the range rhi = [m so \, m atm ] becomes just slightly 
more relaxed. Within a thermal picture of leptogenesis, as we are assuming, the lower 
bound on Mi is closely related to a lower bound on the reheating temperature of the early 
Universe T re h that is approximately given by the expression [7] 

M min 
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Figure 9: Lower bound (3<r) on T re h in the MSSM case (cf. (68)). The lines correspond to 
the same cases of fig 7. For values K > 5 (right of the dashed line) the lower bound does not 
depend on the initial conditions. 

plotted in Fig. 9 for a = 0.5 8 . For K ~ 5 one has the lowest T re h value independent on 
the choice of initial conditions and conservatively we obtain T reh > 10 9 GeV. 

The existence of such a lower bound on T re h represents a problem within the MSSM 
picture, since this has to be compared with the upper bound from the avoidance of the 
gravitino problem. Values of the gravitino mass <C 100 TeV require reheating temperature 
T rc h <C 10 10 GeV [35]. From Fig. 9 one can see that there is compatibility only for values 
mi < m atm . 

These results are valid under the conservative assumption of maximal phase. In gen- 

8 The result is more conservative compared both to the numerical results presented [33] and to the 
analytical ones presented in [32] : the wash-out from scatterings is completely neglected, the approximation 
(62) slightly overestimates (within 10%) the numerical results and the approximation (68) underestimates 
T re h at K < 4. In this way the restrictive bounds we obtain for models with phase suppression, can be 
regarded as very conservative ones. Note also that we are using a different value of a compared to [7] as 
a consequence of the different expression employed for zb{K). 
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eral, if one modifies the picture such that the final predicted asymmetry t\b — ► £,Vb, then 
the lower bound on M 1 and on the reheating temperature change as (M™ in , T r ™ n ) — > 
(M™ m , r™™)/^ [9]- Therefore, if we consider models with some phase suppression, 
sin^ < 1, the minimum allowed values of M\ and T re h increase compared to the case of 
maximal phase. On the other hand, if M™ m becomes larger than about 10 14 GeV, the 
assumption that «f depends only on rh\ is not valid any more, since AL = 2 processes 
give an additional contribution to the wash-out depending also on M 1 [8]. In this case 
the efficiency factor is exponentially suppressed at large M 1 values such that 

~ ~ -— ,yi 

Kf(mi) -> Kf(mi) e z s 1()14 Gev ; (cj ~ 5) (69) 

and ^ ax (Mi, mi = 0, mi) (cf. (34)) has a maximum at Mi ~ 2 z B 10 13 GeV given by [31] 

in-3 e v 

^ ax ( mi = 0, mi) ~ 1.5 x 10~ 5 — (mi > 1(T 3 eV) . (70) 

This implies that the lower bound on the phase (cf. (59)) cannot be arbitrarily evaded 
increasing Mi and one can place an interesting lower bound, independent on the value of 
Mi, given by 

sin 5 L > 4 x 1(T 5 (mi > lO" 3 eV) . (71) 

1U 6 ev 

Let us study two special cases. For p 2 = one has rh\ = ps m atm and from the Eq. (60) 

sin Sl = sin ip 3 . (72) 

One can set 5l = y?3, fixing unambiguously 5l- Again one can see that the effective 
leptogenesis phase is maximal for p 2 = (i.e. Q 2 i = 0) and X 3 = 0, corresponding to 
if 3 = 7r/2. For random assignations of if 3, the average value of the phase, (| sin<p 3 |) = 
2/n ~ 0.6, suggests that a value close to the maximal one is a reasonable possibility not 
requiring any particular tuning. 

Within models with single or sequential RH neutrino dominance, characterized by 
\Xj\ ~ 1 3> l-Xj^jl, there are two choices: either X± ~ 1 3> \X$\ or conversely X 3 ~ 1 > 
|Xl|. The first one corresponds to maximal phase sin^L ~ 1, while m-i ~ i3m atm . If 
furthermore one has I3 < 1, then one has the kind of models we previously discussed, 
with Q obtained as a small complex perturbation of the matrix (25) and small values of 
rhi <^ m a t m , such that minimum values Mi,T re h ~ 10 9 GeV are possible. In the second 
case, i.e. X 3 ~ 1 ^> |-X"i|, such that JVi contributes only to m 3 , one has necessarily some 
phase suppression and moreover rh\ > m atm . The value of the phase can be expressed 
as a function of rh\ and is given by sin 5l = a/1 — m 2 &tni /m\. Since the efficiency factor 
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is approximately inversely proportional to mi (cf. (65)), then one easily finds that the 



M\ and on T re h that can be read from Fig.'s 8 and 9 obtaining M\ > 1.5 x 10 11 GeV 
and T rch > 2 x 10 10 GeV. Moreover the lower bound on the phase (cf. (71)) becomes 
sin 03 > 2 x 1(T 3 . 

Still another possibility is that N\ contributes both to mi and to m 3 and in this case 
l^i I ~ l^sl ^ 1, such that rhi > m atm , still implying M x > 10 11 GeV and T rch > 10 10 GeV. 

The second special case is for p 3 = 0, such that rhi = m 2 p 2 and from the Eq. (60) 
one obtains 



If one considers the case of normal hierarchy, then a = m so i/m atm ~ 1/6 and this un- 
avoidably leads to a strong phase suppression, sinS^ ~ sin<p 2 /6 <C 1, that results in a 
more stringent bound on Mi (see the dashed line in Fig. 8) and on T re h (see the dashed 
line in Fig. 9). For models with N x contributing only to m 2 , such that X 2 — 1 ^> \Xi\, 
the phase suppression is sin^ < 1/6, and similar lower bounds for Mi and T rc h, as in 
the case of Ni contributing only to m 3 , hold, even though the wash-out is much weaker 
because now rhi ~ m so i <C m atm . 

We have thus two extreme cases, p 2 — or p 3 = 0, where the effective phase can 
be maximal or remarkably suppressed. Let us now analyze the general case, with both 
non- vanishing p 2 and p 3 . Since we are interested in the lower bound on Mi, we can study 
the expression (60) for maximal phases, sin(p 2 = sin<p 3 = 1, and thus 



These are the most general models with Ni contributing only to mi, such that \X 2 \, \X 3 \ 
Xi ~ 1. For p 3 ^> crp 2 one has sin 5l — 1 and rhi ~ p 3 m atm , such that the limit p 2 = is 
recovered, while conversely, for p 3 <^ a 2 p 2 , one has the phase suppression sui<5l = a ~ 1/6 
and rhi ~ p 2 m sol and one recovers the limit p 3 = 0. 

In the intermediate case, a 2 p 2 < p 3 < ap 2 , one has 2 a < sin 5l < 1/2. For example 
for p 3 = ap 2 one has rhi = 2m sol ~ 0.016 eV and sm.8 L = 1/2 corresponding to Mi > 
4 x 10 10 GeV and T rch > 7 x 10 9 GeV. For p 3 = a 2 p 2 < 1 one has sin<5 L = 2a ~ 1/3 
and rhi — "^soiP2- The corresponding lower bounds on Mi and T re h are shown in Fig. 8 
and 9 (dot-dashed lines). If p 2 = 1 one gets M x > 3 x 10 10 GeV and T rch > 5 x 10 9 GeV. 
Therefore, this case represents the best compromise to get the lowest possible values for 
and Mi without having p 2 ,p 3 < 1 and pi ~ 1. 




sin 5l = o~ sin <p 2 . 



(73) 



sin<5i = 



P3 + Q- 2 P2 

P3 + 0~p 2 



(74) 
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We can conclude that the possibility to access the window 10 9 < Mi (T re h)/GeV < 
(0.5) 3 x 10 10 relies on a particular class of models, such that p 2 , p 3 <C pi — 1, the only 
ones allowing both maximal phase and rhi <C m atm - This conclusion does not change 
if one considers the case of inverted hierarchy: it is true that one does not have phase 
suppression but values of rhi <C lCT 2 eV are obtained again only for p 2 , p^ <C 1 — p\. 

More generally values T reil < 10 10 GeV are obtained only within models where Aq 
contributes only to m\. Among sequential dominated models these correspond to models 
where m\ is dominantly determined just by Aq. For models where Aq contributes only 
to m 3 , including sequential dominated models with m 3 dominated by Aq, we have seen 
that, if p 2 = 0, then necessarily the lower bounds are more restrictive, such that Mi > 
1.5 x 10 11 GeV and T re h > 2 x 10 10 GeV. Allowing a non zero Y 2 the bounds can only 
become more restrictive. On the other hand if one considers models where Aq contributes 
dominantly only to m 2 (X 2 ~ 1 ^> X lt X 3 ), including those where m 2 is dominated by 
Aq, the limits obtained for p 3 = can be relaxed if one allows a non zero Y" 3 . The 
best optimization between high phase and low rhi is obtained for configurations with 
Y 2 = 0. In this case one has sin^ = 1 — m so \/fhi and it is simply to see that the 
asymmetry is maximized for rhi ~ 2m so i corresponding to sin^L ~ 0.5. One then obtains 
Mi > 4 x 10 10 GeV (within SM) and T rch > 7 x 10 9 GeV (within MSSM) 9 . 

5.2 Quasi-degenerate neutrinos 

In this case one has mi ~ m 2 ~ m 3 and thus mi /rhi — where remember that 
i = pi + p 2 + p 3 > I. The expression (46) can then be re-casted as / ~ \/l 2 — 1 /£ and 

9 An interesting exercise is also to consider the limit Af 3 — > oo, as already studied in [36]. In this 
situation the Q matrix is simply given by 



This implies that Y<i = — 13 and from the eq. (60) one can see that for inverted hierarchy there is a large 
phase suppression. For normal hierarchy, neglecting 0{<j 2 ) terms, one has that sinc>£ ~ fiatm^s/'Wi- 
The maximum of I3 can be found analogously to the general case simply replacing m\ with TO2 in the cq. 
(51), such that Y™ — /(to 2 , mi)/(m atm (1 + a)) (where now rh\ > m 2 ). It is then simple to find that t]b 
is maximized for fh\ ~ 2to so i and that sin<5i ~ 0.4, such that Mi >5x 10 10 GeV (within the SM) and 
T rc h 9 x 10 9 GeV (within the MSSM). These results are more relaxed than in [36], simply because we 
are taking into account the recent results of [33, 7] on the efficiency factor and on the relation between 
M™ ln and T™^ 11 , as previously discussed. 



n = 




(75) 
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X 3 =1/2 




the Eq. (58) for the phase becomes 

2 

sin 5 L ~ — [p 3 sin y? 3 + <x p 2 sin <p 2 ] ■ (76) 

Vt — 1 

For p 2 = and X 3 = 1/2 one has Y3 = 2 / v 7 ^ 2 — 1 and the case of maximal phase is 
recovered. On the other hand if one considers models with p 3 = 0, then the maximum 
asymmetry is still realized for X 2 — 1/2 but this time one has sin 5l = <? 2 , meaning that if 
a = m atm /m so i ~ 1/6 ('normal' quasi-degenerate spectrum), then one has a strong phase 
suppression, a times stronger than in the case of fully hierarchical neutrinos. 

We can also consider an intermediate case with both p 2 and p 3 7^ 0, for example when 
(p 2 = <^3 and X = 1/2 (dot-dashed arrows in Fig. 10). In this case the general expression 
(76) gets specialized into sin 5l = 1 — C + ° r2 C> where ( = p 2 / (p 2 + p 3 ), and describes the 
transition between the case with p 2 = and the case with p 3 = 0. 

Another example is given by a situation where p 2 = but the condition X 3 = 1/2 is 
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not realized. For example if one takes A 3 = (see solid arrows in Fig. 10), then one gets 
sin 5/, = yjl — l/£ 2 . For i >> 1, corresponding to rh\ ^> mi, one recovers maximal phase, 
while for finite values of rh\ there is some phase suppression. An interesting example is 
to take rh\ = \/2mi, since this corresponds to the peak value that saturates the upper 
bound on neutrino masses [7]. In this particular case one finds sin 5l = a/1/2 ~ 0.7. 
An analogous result holds if one takes X± — and A3 = 1. Notice that these two cases 
correspond to a lightest neutrino mass dominated by the lightest RH neutrino (for A 3 = 0) 
or by the heaviest (for X 1 — 0). 

We have also considered the case of configurations with Y2 — and A3 = (1 — A 2 )/2, 
such that when A2 = one recovers the configuration with maximal phase (dashed arrows 
in Fig. 10). In this case one obtains easily 



sin5 L = v /l-^ (77) 

and for example, for p 2 = £/3, one gets siuSl ~ 0.7. These examples show that typically 
a phase suppression is expected. This is relevant in connection with the neutrino mass 
upper bound, m, < 0.1 eV, because it implies that the saturation of the bound actually 
occurs, not only for special values of M 1 and m-i, but also for special QjiS configurations. 



6 A new scenario of thermal leptogenesis 

So far we have assumed that the final asymmetry is not influenced by the two heavier RH 
neutrinos and this has considerably reduced the number of seesaw parameters on which the 
final asymmetry depends. This assumption holds if the asymmetry generated by the A 2 's 
and by the A 3 's is either negligible or is efficiently washed-out by the inverse decays of the 
Ax's and also if the approximated expression (30) for the CP asymmetry e 1 can be used. 
As we will discuss in Section 7, these two conditions require some degree of hierarchy in the 
heavy neutrino spectrum. Here, even though we still assume a hierarchical heavy neutrino 
spectrum, we study a scenario where the final asymmetry is produced by the A 2 's decays. 
This is possible if the value of the effective neutrino mass lies in the weak wash-out regime. 
In this case the wash-out from N\ inverse decays and scatterings is not sufficient to wash- 
out a previously generated asymmetry, in particular that one generated from the A 2 's 
decays. As we discussed in Section 2, the possibility that rh\ < m* ~ 10~ 3 eV, relies on 
models with an Q matrix obtained as a small perturbation of a 23-rotation in the complex 
plane (cf.(25)). Naively one could think that, since we are still assuming a hierarchical 
spectrum of neutrinos such that M 2 ^> M±, and because of the lower bound on Mi, one 
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has to require large values of the reheating temperature such that T rch > M 2 3> 10 9 GeV. 
However, it is easy to understand that actually the lower bound on M 1 disappears and is 
simply replaced by an analogous lower bound on M 2 . First of all let us notice that in a 
general case the final asymmetry is the sum of the contributions from all the three iVj's 
and, assuming a sufficiently hierarchical heavy neutrino spectrum, it can be written as 

T i 

that generalizes the expression (17) with Kf t i = Kf. Let us focus on the asymmetry 
produced by the N 2 decays. These will be described by a quantity, m 2 , that plays an 
analogous role of rh\ in the description of the N\ decays and is defined as 



(m ] D m D 

m 2 = - 



22 



^2 mi ■ ( 79 ) 



M 2 

For our purposes the fl matrix can be conveniently parameterized in the following way 



^ a/1 — ^21 ~~ ^31 ^12 _ \/^21 + ^31 ~~ ^12 ^ 



n = 



1^21 ^22 — \/ 1 — ^22 — ^ 



(80) 



\ ^31 a/1 — fi 22 — Slf 2 V^22 + ^12 ~~ ^31 ) / 



The orthogonality also implies that each product of different columns or rows has to 
vanish. Therefore, we can still express one of the four left complex elements as a function 
of the other three. For example, if f2 21 ^ 1, one can write 10 



fil2 = • (83) 

Since we are assuming pi ~ 1 and p 2 ,p3 <C 1, this expression simplifies into 

Q 12 ~ -n 21 fi 22 ± Q 31 yjl-n"^ , (84) 
10 Alternatively, if f^i = 1 but f^i 7^ 0, one has 

" 12 - 277w^ ' (81) 

In the particular case that fi 2 i = 1 and ^ 31 = 0, then Q 2 2 = and the £1 matrix specializes into 



f2 2 i -y/l - O 



0= I 1 I . (82) 

y/i - n 2 21 n 2 i 
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from which one can deduce that pi 2 < (P21 + P31) (1 + P22), implying that it is impossible 
to have p 12 ~ 1 and p 2 2 <C 1- Therefore, if rh\ <C to*, from the orthogonality of f2 it 
follows that is impossible to have m 2 <C to* too. This means that if the asymmetry from 
the iVi's is generated in the weak wash-out regime, then the asymmetry from the iV 2 's has 
necessarily to be generated in the strong wash-out regime. This is an interesting result 
because, if the final asymmetry is generated by the iV 2 's, then it will be independent on 
the initial conditions exactly like in the typical case when the final asymmetry is generated 
by the iVi's and rh\ lies in the strong wash-out regime 11 . 

Moreover notice that a possible asymmetry generated by the heaviest RH neutrinos, 
the iV 3 's, is always washed-out and negligible within the assumption of a hierarchical RH 
neutrino spectrum. 

We have now to calculate the final asymmetry generated by the iV 2 's and see whether 
this can explain the observed baryon asymmetry. Therefore, we have to calculate the 
CP asymmetry associated to the N 2 decays. The expression for e 2 is simply given by an 
expression similar to the Eq. (26), 



Im [(fet h )l] 

£ 2 - -7T- > , 7TTT\ ~ x 



8vr E {Wh) 22 



This time however it cannot be written in the form (30), since now the two terms have 
to be evaluated in two different limits: the first one (i = 1) for Mf/M 2 <C 1 such that 



f ^Ml) +fs {Ml)^'M- 2 { ln M- 1 - 2] (86) 



and the second one (i = 3) in the limit .Ut,' .1 /.; 3> 1 such that 



, Mo 2 \ , /Mo 2 \ 3 M 2 



In these limits the Eq. (85) becomes 

~3M 2 r/1 +IN2 . Mi /, M 2 



e 2 ~ 



8tt (/it ^) 22 



2M, 



M(rfh)l 3 ] + ^ (in ^-2) Im[(/|t/i) 



(88) 



11 This result is also interesting because it suggests that after all an initial asymmetry is always washed 
out even if the lightest RH neutrino decays occur in the weak wash-out regime and so that from this 
point of view there is no problem with the initial conditions. However, in this case one has to impose 
that the reheating temperature is higher than ~ M2I zb(K2) and so the lower bound on T rc h would be 
much more restrictive. Moreover the problem of the dependence on the initial N\ abundance in the weak 
wash-out regime still remains. 
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From the Eq. (5) one can easily derive the relation 



(/it h) 



Yl m hW hi n hj , 



(89) 



that can be used to recast the Eq. (88) in terms of the f2 matrix elements and of m 2 as 
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(90) 

Let us now consider three different limit cases of models where pi — 1 and p 2 , p 3 <C 1. The 
first case is for p 2 = and p 22 = 1, corresponding to a small complex angle 13- rotation 



\ 





^31 



-Q31 \ 

1 

xA 3 ^/ 



(91) 



As we know, in this case the CP asymmetry E\ is maximal. On the other hand it is 
simple to see from the Eq. (90) that £2 = and thus the asymmetry from the AVs 
vanishes and the traditional picture still holds. The second case is for p 3 = and p 22 — 1, 
corresponding to a small complex angle 12-rotation 



( - n 2 21 -n 21 

Q 21 y f l^W l 
\ 1 



(92) 



From the Eq. (90) one can easily see that the first term vanishes and one has 



£2 



Mi m atm M 1 
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a ( In ^ - 2 ) sin<p 2 = 77 £ m ax(Mi) a ^ ( In ^ - 2 ) sin<p 2 . (93) 



Mo 



Mo 



Mi 



This result has to be compared with e\ = £ max (Mi) a sin<p 2 (cf. (31), (32) and (60)) 
showing that e 2 < £i- Therefore, in this second case too, one has that the asymmetry 
generated from the iV 2 's can be neglected and again the usual picture still holds. 

The most interesting case is for p 2 = p 3 = 0, corresponding to a complex 23-rotation 
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(94) 
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The second term in the Eq. (90) vanishes while the first, after some algebraic manipula- 
tions 12 , gives 

3 M 2 ml- ml lm(Ql 2 ) 
S2 ^ — I , 95 

16 7T IT m 2 

showing that now e 2 can be large, while on the other hand, as we know, e± — 0. Therefore 
we see that, when rh\ <C m±, the two possibilities of having maximum CP asymmetry for 
the Ni decays or for the N 2 decays are complementary and when one is maximum the 
other vanishes and vice versa. 

The maximum value of \e 2 \ in the Eq. (95) can be found analogously to when we 
maximized S\ for configurations with Q 2 i = 0. Now the role of and flfi is replaced by 
Vt 22 and Vl\ 2 and one has to find the maximum of Y 22 = Im(fi| 2 ). This can be calculated 
solving a set of two algebraic equations for F 2 m and for the corresponding value X 22 = 
Re(f2 2 2 ), exactly equal to the Eq.'s (49) and (50) but with the obvious replacements 
m\ — > m 2 and rh\ — > rh 2 . We can also introduce again a function f(m 2 ,rh 2 ) related to 
Y 2 ™ exactly as the function /(mi, mi) was related to Y m (cf. 51), writing 

t( ~ ^ m 2 + m 3 

f{m 2 ,m 2 ) = Y 22 < 1 (96) 

m 2 

and in this way one can write for the maximum CP asymmetry 

eT x = £max(M 2 ) ms ~ m2 f(m 2 , m 2 ) . (97) 

m atm 

Notice that this is just the bound (43) for e±, where now m± is replaced by m 2 that is fixed, 
except for the possibility to choose between normal and inverted hierarchy. Therefore, this 
time e™ ax depends on whether one considers normal or inverted hierarchy. The maximum 
value is obtained for normal hierarchy because in this case m 2 m 3 and also because the 
function f(m 2 , rh 2 ) cuts off the asymmetry for rh 2 < m 2 , that means for rh 2 < m so \ in the 
normal case and for rh 2 < m atm in the inverted case. One can understand this difference 
also considering that in the case of inverted hierarchy one has a degeneracy m 2 ~ m 3 and 
thus it is analogous to the case of a quasi-degenerate spectrum for E\. 

The bound (97) on e 2 will give rise to a lower bound on M 2 analogously as the bound 
on E\ determined a lower bound on Mi (cf. (61)). The difference is that now the function 
f(m 2 ,rh 2 ) is less than one, even though we are dealing with hierarchical neutrinos and 
thus one has for the SM (MSSM) case 

M 2 > y)xlO»_GeV 
f{m 2 ,m 2 ) Kf{m 2 ) 

This lower bound on M 2 is shown in Fig. 12 for both normal (thick solid line) and 



12 In particular we used that 2 Re(fi2 2 \A — ^22) im (^22 \/l ~ ^22) = ~ I m (^22)- 
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Figure 11: lower bounds (3<r) on M 2 (solid lines) and on r re h (thin lines) obtained within 
models with rhi <C m+ and P2 = P3 = 0. The solid lines refer to the case of normal hierarchy 
while the dashed lines refer to the case of inverted hierarchy, (left) SM case; (right) MSSM case. 

inverted (thick dashed line) hierarchy. Obviously, there will be again an associated lower 
bound on the reheating temperature given by the Eq. (68), where now M™ m has to be 
replaced by M™ in . This lower bound on T re h is also shown in Fig. 11 (thin lines). The 
lower bounds on M2 and on T re h can be calculated again both within the SM (Fig. 11 
left) and the MSSM framework (Fig. 11 right). 

Notice that within this scenario, with P2 = P3 = 0, there is no lower bound on the 
value of Mi, that thus can be arbitrarily small. This occurs because we are assuming 
rh\ < m*, while if rhi > 5 m* the asymmetry from the iV 2 's would be fully washed-out 
from the N\ inverse decays and in this case there would not be any allowed value for 
Mi. Therefore, for P2 = P3 = 0, there is a transition between a regime where the lower 
bound vanishes, for rh\jm^ -C 1, and a regime where the lower bound goes to infinite, for 
rhi > 5 m*. In the region < rh\ < 5 m* there is quite a sharp transition that can be 
determined by a detailed description of the wash-out of the asymmetry produced by the 
AVs from N\ inverse decays. We have not calculated the lower bound in this transition 
region, that in Fig. 8 is schematically indicated with a box. 

Since we are assuming P3,p2 «C p\ — 1, then approximately m 1 ~ v 2 \h u \ 2 /M 1 . 
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Therefore, if Mi becomes as low as the electroweak scale and because mi <C lCT 3 eV, 
then one has \h n \ <C 1CT 7 . This corresponds to a situation where the lightest right 
handed neutrino is decoupled and its interactions are switched-off, with no possibility to 
detect it in the accelerators 13 . 

(2) 

Notice that one could also introduce an effective leptogenesis phase, sin 8 L , associated 
to the asymmetry generation from the A^'s and study its suppression when one considers 
cases different from the maximal production. For example this happens if Y22 7^ Y m , 
but also if one considers a combination of the three limit cases that we have studied, 
described by an Q matrix that is the product of the three complex rotations. In such 
a general situation one can have both a production from the Nis in the weak wash-out 
regime and a production from the AVs in the strong wash-out regime. However, not to 
have again problems with the dependence on the initial conditions, one should impose 
that the contribution from the A^'s is the dominant one, even though not maximal. 

To summarize, we have seen that the possibility of an asymmetry generated by the A^'s 
is realized for quite a special form of the matrix Q (cf. (94)). Nevertheless this scenario is 
interesting because it solves the problems associated with the generation of the asymmetry 
from the Aq's with rh\ in the weak wash-out regime. Indeed, we have already shown that, 
since m 2 > m sol ^> m*, there is no dependence on the initial conditions, similarly to the 
typical case of asymmetry generated by the iVi's with rh\ > 5 m* 14 . Another interesting 
point is that now the value of rh\ has not to be fine-tuned such that mi < mi <C m* 
but one has simply rh\ — mi. Note that the leptogenesis conspiracy [31] applies also to 
this scenario. This is the observation that the measured atmospheric and solar neutrino 
mass scales lie in the correct range of values, between 10~ 3 eV and 1 eV, for leptogenesis 
to work. Indeed if m so i and m atm were found much larger than 1 eV, then m 2 > m so i 
would have yielded a too strong wash-out. On the other hand if they were found much 
smaller than 10~ 3 eV, then the maximum CP asymmetry would have been equally smaller 
and the lower bound on T rch at least two orders of magnitude more restrictive. Moreover 
the asymmetry production would have occurred in the weak wash-out regime with the 

13 Unless one introduces some extra gauge interaction, as recently proposed in [37] within the context 
of resonant leptogenesis. In our case it would be even simpler, since no ad hoc motivation for a strong 
degenerate heavy neutrino spectrum is needed. 

14 Notice that we could have also considered another possibility with mi in the strong wash-out regime 
and ffi2 in the weak wash-out regime such that the N 2 decays occur below the temperature Tg where the 
wash-out processes involving the AVs freeze out. In this case the wash-out would be circumvented and 
the final asymmetry can also be generated by the AVs decays or even by the A^'s if fri3 is in the weak 
wash-out regime instead of fh 2 - However, this possibility simply transfers the problem of the dependence 
on the initial conditions from the AVs to the AVs (or to the AVs). 
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related problems of dependence on the initial conditions. In conclusion, even though 
this scenario sits in some special corner of the seesaw parameter space, it represents an 
interesting possibility to be taken into account. 

7 On the hierarchy of the heavy neutrino spectrum 

Our results have been obtained under the assumption that the spectrum of the heavy 
neutrinos is hierarchical and this justified some approximations in the calculation of the 
final asymmetry. We want now to discuss the conditions of validity of these approxima- 
tions. A first aspect concerns the possibility to neglect the asymmetry and the wash-out 
from the decays and inverse decays of the two heavier neutrinos. This is possible if M 2 is 
large enough compared to Mi. A conservative assumption is to impose that the decays 
and inverse decays of the two lightest RH neutrinos do not interfere with each other such 
that the generation of the asymmetry from the iVi's and that from the iV 2 's proceed inde- 
pendently. This holds if the generation of the asymmetry and the wash-out from decays 
and inverse decays of the iVi's start only after the end of the analogous processes from 
the 7V 2 's. 

Under this assumption and if fhi lies in the strong wash-out regime, the rate of gen- 
eration of the asymmetry from the iVi's, dr]B/dz, is well approximated by a Gaussian 
centered around a peak value zb(Ki) [5, 7], where Ki = rhi/m*. More explicitly the 
efficiency factor can be written as a Laplace integral given by 

«£,!=/ dz'e-^ z,) ~ dz'e-^ z '\ (99) 

where zb(Ki) is that value of z' where the function ip(z'), determined both by the decay 
rate and by the wash-out from inverse decays, has a minimum. A second order expansion 
around zb(Ki) makes possible to approximate the exponential with a Gaussian centered 
at zb(Ki). This procedure is valid only if the inverse decays from the iV 2 's are inefficient 
for z > zb(Ki) — A(Ki), otherwise there would be an additional contribution to /t f l 
coming from the N 2 inverse decays wash-out. If m 2 lies in the strong wash-out regime 
too and if one neglects the wash-out from JVi inverse decays 15 , then the generation 
of the asymmetry from the iV 2 's will be also described by a Gaussian centered around 
a value z ~ z B (K 2 ) Mi/M 2 , with K 2 = m 2 /r?v This means that for values of z ~ 
[z B (K 2 ) + A(K 2 )] Mi/M 2 the generation of the asymmetry from the N 2 decays and the 

15 This is a conservative assumption, an account of this effect would go into the direction of relaxing 
the final condition on the ratio M 2 /Mi. 
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wash-out from the N 2 inverse decays can be neglected. Therefore, they will not influence 
the final value of the asymmetry, if the peaks of the two Gaussian are sufficiently separated, 
that means if [z B (K 2 ) + A(K 2 )} M 1 /M 2 < z B {K x ) - A(Ki). 

Let us impose a separation of the two peaks of about 3 a that guarantees a precision 
much less than 10% in the calculation of the asymmetry when N 2 decays and inverse 
decays are neglected. Using the expression (64) for zb{K) and considering that the width 
of the Gaussian is approximately constant with K and given by a ~ 1.5 [7], this implies to 
have A(Ki) ~ A(K 2 ) ~ 1.5 a ~ 2. Therefore, we arrive to the following very conservative 
condition 

jft + 2 . (100) 

Mi " ZB^Kr) - 2 y ' 

Taking the most conservative choice of values 16 , K x ~ 10 and K 2 ~ 10 4 , the RH side 
of this inequality is maximized and one finds that: if M 2 /M\ > 5, then the effect of N 2 
decays and inverse decays on the final asymmetry can be safely neglected. If fhi is in the 
weak wash-out regime, then one can have a situation as described in Section 6, while if 
m 2 is in the weak wash-out one can have a situation as described in footnote 14. 

The second aspect to be considered is the CP asymmetry and in particular the condi- 
tions of validity of the approximations that from Eq. (26) lead to Eq. (30). The first one 
can be written as 

where Xi = Mj /Ml and the function £(x), given by [9] 

1 + x\ 2 — x 



2 

= -x 



(1 +x) In 



(102) 



x ) 1 — x 

approaches 1 for x ^> 1. After some easy algebra, the Eq. (101) can be re-casted as 

ei = £(x 2 ) ei(Mi, mi, mi, Q 2 n ) + [£(x 3 ) - £{x 2 )] Aei , (103) 

showing that in the limit x 2 , x^ — > oo one recovers the result S\ ~ £i(Mi, mi, mi, Q^) (cf. 
30), while, for finite x 2 , x%, there is both an enhancement £,(x 2 ) of the usual term, plus an 
additional contribution [£(x 3 ) — £,(x 2 )] Ae±, where 

_ 3 Im[(fet h)j 3 ] 1 MEa ^^d 2 nn4 , 

16 71 {hm) u yfxz m atm mi 



16 This choice corresponds to have N\ decays in the strong wash-out regime and rri2 ~ 10 eV. This 
large mi value can be obtained if mi ~ m a tm ^32! an d ^32! ~ 10 3 - Since M3 > 10 10 GeV, from the eq. 
(6) one can check that values rh^ S> 10 cV would imply Yukawa couplings much larger than 1. 
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In the second expression we have used the relation that connects h to Q (cf. (89)) and 
the definition of e max (cf. (41)). In this way we can write 



£ £ = ~7~/T/FT = eWWrm,^) sin5 L + K(* s ) - *(* 2 )] ^gi^] _ (1Q5) 



First of all let us notice that in this general case the asymmetry depends on 4 additional 
parameters, the masses M 2 and M 3 and the two Q parameters that were disappearing 
within the assumption of heavy hierarchical spectrum. If we use the parametrization Eq. 
(80) for Q, then these two additional parameters can be identified with the real and the 
imaginary part of Q 22 . 

It is instructive to study £ £ in two limit cases. For x 2 = ^3 (i-e. M 2 = M 3 ) one 
has Aei = and the enhancement of the asymmetry, compared to the hierarchical case, 
is simply described in terms of £(x 2 ) [9]. Therefore we see that in this case the final 
asymmetry will depend just on one additional (seventh) parameter compared to the case 
when a full hierarchical heavy spectrum is assumed. Like for the effect of asymmetry 
and wash-out from the iVVs, one has again that if M 2 > 5 Mi, then the precision of the 
approximation of heavy hierarchical spectrum is much below 17 10%. 

The second limit case is obtained for x 3 — > 00 such that the difference £(x 3 ) — £(x 2 ) c± 
1 — £(^2) is maximum with respect to x 3 . It is useful to calculate Ae 1 in the case of 
maximal effective phase and full hierarchical neutrinos such that Q21 = X 3 = 0; the 
general expression (105) becomes simply 



One can see that a crucial role is clearly played by the imaginary and by the real part of 
f2 22 . If X 2 2 + Y 3 Y 22 = —1, one has the curious result that S\ = e max (Mi), irrespective of 
the value of x 2 . Thus we see that the enhancement of the CP asymmetry, when M 2 — > M 1; 
is not unavoidable, but a possibility. In general, if X 22 + F 3 Y 22 ~ 0(1) and if M 2 > 2 M 1 , 
one has that the first term in the (101) dominates, while if X 22 + Y 2 Y 22 is much larger 
than one, then a dominance of the second term is possible. However, if one imposes 
M 2 > 5 Mi, large values £ e > 2, necessary to evade significantly the lower bounds on Mi 
and on T reh , rely on unnaturally huge value of X 22 + Y 3 Y 22 , implying huge values of rhi 
or of m 2 or both. Analogous results have been obtained also in [38] . 

Notice that in general one can say that the two particular cases we have studied, 
x 2 = x 3 and x 3 = 00, can be regarded as limit cases of two different classes of RH 
neutrino spectrum: inverted heavy neutrino spectra, with Mf — M| < M| — M\ , and 
17 For M 2 =5Mi, one has £(x 2 ) - 1 ~ 0.02. 
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S e = Z(x 2 ) + [S(x 2 )-l][X 22 + Y 3 Y 22 ]. 



(106) 
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normal heavy neutrino spectra, with M| — Mf > M| — M\. For small imaginary parts 
and small values of Ae x the second term in the Eq. (105) is small and £ e is the same 
in both cases, while for large imaginary parts the case of normal heavy neutrino spectra 
gives rise to a larger enhancement arising from Ae±. However as far as M 2 > 5 Mi the 
enhancement can be relevant in changing considerably the lower bounds on Mi and on 
T rc h only for large values of the imaginary parts. 

Another issue concerns the effective leptogenesis phase. The second term in the Eq. 
(105) is not maximal for the same configurations that maximize the first one and thus the 
effective leptogenesis phase sin Sl is related only to the first term. In particular we have 
seen that the second term is sensitive also to the values of X 22 and Y 22 . Here we are not 
interested in a systematic analysis of Ae±, but it is interesting to describe two particular 
cases. First of all notice that if rhi = mi, that is equivalent to assume an Q matrix as in 
the Eq. (25), then not only the first term vanishes, because /(mi, Si) = 0, but one can 
check that Aei = too. This is important because, otherwise the results that are usually 
presented in the limit fhi — > mi, would have been valid only under special conditions. On 
the other hand we can consider a simple case showing that when sin^L = the second 
term does not necessarily vanish. Let us take X 2 — Y 2 — Y 3 — 0, such that sin5 L = 0, 
and let us also assume for simplicity the case of full hierarchical neutrinos, such that 
rhi = m 3 X 3 . In this case one can easily arrive to the expression 

& = [tM-tM](l-X 3 )Y 22 , (107) 

showing that the second term does not vanish in general when sin Sl = 0. However, once 
more, we can see that if M 2 > 5 Mi, then one has to require Y 22 ^> 1 to have £ £ > 1. 

Summarizing, we can say that if the condition M 2 > 5 Mi holds, then the usual 
approximations work very well with a great precision, except for some situations involving 
large imaginary parts of the fi^'s. 

If one considers quasi-degenerate light neutrinos with mi > m atm , then, while the first 
term in the (105) gets suppressed as oc (m 3 -mi)/m atm , the second does not and therefore 
it can become dominant more easily, especially if one considers a 'normal' heavy neutrino 
spectrum with x 3 ^> 1 and finite x 2 [19]. This is relevant in connection with the neutrino 
mass upper bound. One has to say however that the possibility of having mi > m atm holds 
only for a restrictive choice of the parameters. In particular for reheating temperature 
T rch <C 10 11 GeV the bound falls in the hierarchical regime and thus one does not expect 
important corrections, but further investigation is needed on this point. 
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8 Summary and final discussion 



Leptogenesis is not only an attractive way to explain the baryon asymmetry of the Uni- 
verse but also a powerful cosmological tool to get a unique information on the seesaw 
parameters. It is intriguing that the CP asymmetry underlying the origin of the matter- 
anti matter asymmetry of the Universe could be related to a sort of 'dark side' of the 
seesaw parameter space, the orthogonal matrix Q, that escapes Earth experiments but 
leaves a cosmological imprint. 

Our investigation could cast light into this dark side thanks to a geometrical repre- 
sentation of the orthogonal seesaw matrix, resembling the use of unitarity triangles in 
the study of CP violation in the weak interactions of quarks. Interesting features of 
leptogenesis and of its implications for neutrino mass models have emerged. 

With the use of the seesaw geometry it has been possible to calculate the CP asym- 
metry bound that gives rise to the leptogenesis constraints on the neutrino masses and 
on the reheating temperature, fully characterizing the models that saturate the bound. 
At the same time it has been possible to study the phase suppression in general models 
compared to the case of maximal asymmetry and in this way we could place an interesting 
lower bound (cf. (59)). 

The results have specified which neutrino mass models can access different range of 
values of Mi and T rc h. In particular, 

• the lowest allowed range for Mi,T re h ~ (10 9 — 10 10 ) GeV is actually accessible only 
to a particular class of neutrino mass models where the lightest neutrino mass is 
dominated by the contribution of the lightest RH neutrino, a class of models that 
cannot be easily motivated within models where neutrino masses resemble quark 
mass models [25] 18 ; 

• models with the lightest RH neutrino Aq dominantly contributing to one of the two 

heavier (light) neutrino masses, m 2 and m 3 , that are typically considered to explain 

large mixing angles resembling quark mass models, undergo both phase suppression 

and strong wash-out such that the lower bounds become more restrictive: M\ > 

1.5 x 10 11 GeV (within SM), T reh > 2 x 10 10 GeV (within MSSM) if Aq dominantly 

contributes to m 3 and M 1 > 4 x 10 10 GeV, T rch > 7 x 10 9 GeV (within MSSM) if 

Aq dominantly contributes to m 2 . 

18 This possibility can be however realized in more sophisticated models of fcrmion masses, as for 
example in a recent 6 dimensional orbifold 50(10) model [39]. 
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We have also studied whether these restrictions can be circumvented when the usual as- 
sumption of hierarchical heavy neutrino spectrum is dropped. A significant enhancement 
of the final asymmetry relies on the possibility either of a very fine tuned degeneracy such 
that M2 — Mi, either on large values of the imaginary parts of the Q parameters, we have 
seen the special role played in particular by Q22, or 011 some combination of the two things. 
A better understanding of the relevance of these models is certainly needed, however at 
the moment it seems that their feasibility relies on some ad hoc model engineering. 

In conclusion the results we have obtained can be regarded either as something that 
exacerbates the problems of thermal leptogenesis related to large values of Mi and of 
T re h, the most popular being the gravitino overproduction within traditional supergravity 
models, either as a support of thermal leptogenesis to less typically considered possibilities 
for neutrino models or, a third possibility, as a support to inflationary models with quite 
large reheating temperatures. 

Further investigations are needed to understand which one of these possible interpre- 
tations is the correct one. If it will prove possible to built realistic models with the lightest 
neutrino mass dominated by the lightest RH neutrino, then our analysis has revealed also 
another interesting possibility, namely that the final baryon asymmetry could be actually 
produced not from the decays of the lightest RH neutrinos, the N^s, as usually consid- 
ered, but from those of the second lightest, the iV 2 's. This is a particular scenario working 
for rhi < m+ (cf. (94)), that however does not suffer of the problem of the dependence on 
the initial conditions. 

In conclusion, our analysis confirms that thermal leptogenesis is a very predictive 
model of baryogenesis, implying tight and non trivial correlations between the final baryon 
asymmetry and neutrino mass models. It is impressive the amount of implications that is 
possible to derive by requiring the explanation of just one single observable, the baryon 
asymmetry, within an eighteen parameter model like the seesaw. A large variety of neu- 
trino mass models can explain the baryon asymmetry within the seesaw mechanism but 
non trivial constraints apply and the value of the reheating temperature emerges as one 
of the most crucial parameters on which unfortunately we have little experimental in- 
formation. Future investigations should be hopefully able to understand whether these 
constraints will somehow point to a particular neutrino mass model, able to explain all ob- 
servations simultaneously, from neutrino experiments to the baryon asymmetry, or some 
severe conflict will arise whose solution will require one of the many possible extensions of 
the minimal leptogenesis model with the problem to discriminate in favor of one of them. 
At the moment it is quite remarkable how the oldest minimal leptogenesis model is still 
alive, having successfully passed important tests represented by the values of the neutrino 
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mixing scales, that well allows to talk of a leptogenesis conspiracy [31] and with the future 
perspective that this conspiracy can become even tighter if the absolute neutrino mass 
scale upper bound, rrii < 0.1 eV, will be confirmed by the experiments. 
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Appendix 

Let us show that starting from a generic configuration it is always possible to find, for a 
fixed value of mi, a configuration with higher value of (3 and Y 2 = 0. There are two cases 
to be distinguished. 

The first case occurs starting from a configuration with Y 2 < (i.e. Y < Y 3 ). An 
example of such a configuration is shown in Fig. 3 with dotted arrows 19 . It is then simple 
to construct a configuration with Y 2 = (Y = Y 3 ), same value of rh\ and higher value of 
(3 (or equivalently of the CP asymmetry). Let us fix fi^, that means the values of X 3 and 
Kj, and let us pass to a configuration with a smaller p 2 . In this case one has that the sum 
P2 + Pi = 4/ m i ~ P2 (m 2 /mi — 1), with constant A 3 , where Aj, = rh\ — Piin^, increases. 
This implies that decreasing p 2 with a continuous transformation, one necessarily arrives 
to a configuration with Y 2 = 0. 

The second case occurs starting from a configuration with Y 2 > 0. An example is 
shown in Fig. 3 with dashed line arrows 20 . Let us fix this time f^, that means the 
values of X and Y. In this case one can again pass to a configuration with smaller p 2 , 
same value of rh\ and higher value of /3. However in this case one can end up either again 
with a configuration with Y 2 = 0, shown in Fig. 3 with the solid line arrows, or with a 
configuration such that Q 3l and are aligned, corresponding to ip 3 = ip 2 and shown in 
Fig. 4 with dotted arrows. This because if one decreases p 2 keeping fhi constant, then the 
sum p 2 + p 3 = A\jm 3 + p 2 (1 — m 2 /m 3 ) necessarily decreases while, as we will show in a 
moment, Y3 increases and this implies that j3 increases too. Let us show that Y3 increases 
for decreasing p 2 . It is easy to find an expression of Y 3 as a function of p 2 using simple 

19 Notice that one can always assume \<p2\ < \<P3\, since this has always higher value of (3 compared to 
a dual configuration with \if2\ < \<fi3\ obtained from a rigid rotation of 0|i + ^11 around the axis that 
passes through the points (X 3 ,Y 3 ) and (1,0). 

20 This time, analogously to the first case, one can always assume ^3 > y>2- 
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geometrical relations, obtaining 



Y(X 2 + Y 2 + p Ap) ± X y/[p 2 - (X 2 + Y 2 )] [X 2 + F 2 

2(x 2 + r 2 ) 



(108) 



where we defined p = p 2 + Ps and Ap = p 3 — p 2 . If F 2 < 0, as we are assuming, it is 
then easy to show that dY 3 /dp 2 < — 2F(p — Ap) < 0, showing that Y 3 increases when p 2 
decreases. 

If the initial configuration is such that Ai/m 3 > X 2 + Y 2 , then it is possible to decrease 
Y 2 until this vanishes, otherwise one ends up with a configuration such that <p 2 = ip 3 . In 
this second case one can however perform a different kind of transformation fixing this 
time the values of X and X 3 (other than that of mi) and if one decreases Y 2 then (3 
increases. This can be shown easily taking the derivative of the Eq. (cf. (39)) with 
respect to Y 2 and finding an expression for dY 3 /dY 2 given by 



Thus the CP asymmetry increases when Y 2 decreases. One can try to decrease Y 2 until 
it vanishes. In this case one can again apply the previous result valid for configurations 
with Y 2 = 0. However this is not guaranteed to be possible, since the result (110) is valid 
only for an infinitesimal transformation starting from <p 2 = <p 3 and it is not in general 
possible to decrease Y 2 of a finite arbitrary quantity, while the CP asymmetry increases 
and rh\ is kept constant. However in any case such a transformation brings again to a 
configuration with <p 3 > ip 2 . One can then again shorten p 2 through a transformation with 
fixed values of Y and X and again this can bring either to a configuration with Y 2 = 
or to one with <p 3 = <p 2 . In the second case the procedure can be repeated recursively 
until a configuration with Y 2 = is obtained anyway. In Fig. 4 we give a very schematic 
example of such a procedure. We have thus demonstrated that, starting from any generic 
configuration, there is always another configuration with the same value of fhi, higher 
value of P (or equivalently higher value of the CP asymmetry) and Y 2 = 0. 
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